Propagation of sound in a Bose Einstein condensate in an optical lattice 
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We study the propagation of sound waves in a Bose-Einstein condensate trapped in a one- 
dimensional optical lattice. We find that the velocity of propagation of sound wavepackets decreases 
with increasing optical lattice depth, as predicted by the Bogoliubov theory. The strong interplay 
between nonlinearities and the periodicity of the external potential raise new phenomena which are 
<-h ' not present in the uniform case. Shock waves, for instance, can propagate slower than sound waves, 

due to the negative curvature of the dispersion relation. Moreover, nonlinear corrections to the 
■ Bogoliubov theory appear to be important even with very small density perturbations, inducing a 

^^^^ ' saturation on the amplitude of the sound signal. 
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The study of Bose Einstein condensates in optical lat- 
tices is a very active field of research, both from the the- 
oretical and experimental sides. The presence of the lat- 
tice can drammatically modify the behaviour of the sys- 
tem with respect to the uniform case, giving rise to new 
phenomena like, for instance, a Mott-insulator superfluid 
phase transition hi |2| or the occurence of dynamical in- 
stabilities 0,0>IElE0il3- Several efforts are also focusing 
on the creation of quantum computers ||| and ultrasen- 
sitive interferometers |10). In this paper we study the 
propagation of sound waves on top of a Bose-Einstein 
condensate at rest in a one-dimensional lattice. 

The propagation of sound in a harmonically trapped 
condensate without lattice has been already observed ex- 
perimentally 0, and studied theoretically 0, 0, 
[TEL Il6j . Generally speaking, it is important to study 
the propagation of sound also in the non linear regime, 
where density fluctuations are comparable to the back- 
ground density of the condensate. The reason is that 
only large amplitude wavepackets can be realistically ob- 
served experimentally. It turns out, however, that non- 
linear effects in the sound propagation are particolarly 
interesting also from the theorethical point of view. For 
instance, the formation of shock waves in front of a 
bright sound wavepacket (positive density variation) a 
uniform or harmonically trapped condensates has been 
studied in Hq. and analitycal solutions have also been 
discussed in |l7l Il8l Il9j . In the presence of a peri- 
odic potential, theoretical aspects of the propagation of 
sound have been so far explored o nly in the linear regime 
0, E3, ISH 113, E3, [2^, |23, [23, 113, H3|. The main 
question addressed in this work is the strong interplay 
between nonlinear effects and the periodicity of the ex- 
ternal potential. 

We assume a fully coherent condensate described by 
an effective one-dimensional Gross-Pitacvskii (GP) equa- 
tion. In the uniform case, the spectrum of elementary 
excitations is given by the well-known Bogoliubov disper- 
sion relation u>(q) — yj q 2 /2m (q 2 /2m + 2U), where q is 
the momentum of the excitation, m the atomic mass and 
U = gn the inverse compressibility with g = Airh 2 a/m 
the interaction strength and n the density. For mo- 
menta q smaller than the inverse of the healing length 



£ = 1/ \/8iran, the excitation spectrum depends linearly 
on the momentum with a slope c = yJU/m which defines 
the sound velocity. 

We now consider a condensate in a one dimensional 
optical lattice 



Vop t (x) = s E R sin 2 (j-j-j 



(1) 



characterized by a lattice spacing d and a depth s in units 
of the recoil energy Er = Ti 2 ir 2 /2md. The elementary ex- 
citations of a condensate in presence of a one dimensional 
lattice have been discussed in detail, e.g., in j^- One of 
the main features is the formation of a band structure in 
the Bogoliubov spectrum when the lattice is turned on, 
in analogy with the linear Bloch theory. The spectrum is 
periodic in the quasi-momentum q of the excitations and, 
for a given q, many excitation energies, labeled by a band 
index, are available. Energy gaps open at the boundary 
and at the center of the Brillouin zone, given respectively 
by the Bragg momentum qs = hir/d and q R = 0. 

As in the uniform case, the energy of low energy exci- 
tations is linear in the quasi-momentum, Kuj = c\q\, and 
the system admits sound waves with velocity 



(2) 



The inverse compressibility U increases very weakly with 
the optical lattice depth s, and the effective mass m* in- 
stead strongly increases with s, reflecting the correspond- 
ing quenching of the tunneling rate through the barriers. 
As a consequence the sound velocity is predicted to de- 
crease for increasing lattice depth. In the tight binding 
limit, the effective mass is related to the tunneling pa- 
rameter 8 by to* jm = 2Er/ {it 2 8) and the sound velocity 
can thus be written as c — dy/US/h |2^l2flj. 

However, it is not obvious whether a sound signal of 
observable amplitude can propagate also in deep lattices, 
where the tunneling rate is very small and nonlinear 
corrections to the Bogoliubov theory can be important. 
For instance, nonlinear effects combined with the pres- 
ence of periodic trapping potentials drammatically affect 
the center-of-mass motion of harmonically driven conden- 
sates due to the onset of dynamical instabilities 0, 0, • 
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This paper is organized as follows: In Sect. [I] we out- 
line the procedure to generate a sound signals. In Sect.HTI 
we introduce our theoretical analysis, based on both the 
Gross-Pitaevskii (GP) equation and the discrete nonlin- 
ear Schrodingcr equation (DNLS). In Sect. Illll we obtain 
the velocity of propagation of large amplitude wavepack- T 
ets, and in Sect. llVl we discuss the different regimes for 
sound propagation. In Sect.[V]we comment on the conse- g 
quences of an additional harmonic trapping and, finally, 
in Sect. IVll we conclude discussing the experimental ob- 
servability of the predicted effects. 



I. GENERATION OF SOUND SIGNALS 

One possibility to create excitations in the condensate 
consists in turning on and/or off a perturbation poten- 
tial in the center of the system, for instance by means 
of a far detuned laser beam. If the width of the per- 
turbing potential is much larger than the lattice spac- 
ing, only quasi-momenta much smaller than the Bragg 
momentum are addressed. This procedure generates 
a pair of wavepackets propagating symmetrically out- 
wards. Moreover, only Bogoliubov bands with energy 
lower than the inverse time scale T p of the perturbation 
are excited. In particular, since the gap between first 
and second Bogoliubov band at the center of the Bril- 
louin zone is of the order of 4£ , fl, one has multiple-band 
excitations if T p < K/4Er, otherwise if T p 3> h/iEp 
only the lowest band is addressed. The result of a fast 
perturbation (T p = h/Ep) at low lattice depth is shown 
in Fig^ One finds two pairs of wavepackets propagat- 
ing at two different velocities. The velocity of the slower 
ones is given by the sound velocity @, while the veloc- 
ity of the faster ones is given approximatively by 2qp/rn, 
close to the derivative of the second and third band of 
the spectrum at small q. Those fast wavepackets are not 
sound waves but are superpositions of single particle ex- 
citations: they disappear in absence of the lattice, but if 
the lattice is present they are found in a non-interacting 
gas as well. Apparently, they travel without changing 
shape, simply because the curvature of the spectrum is 
too small to observe their dispersion on the time scale 
of our simulation. This effect can be observed only in 
presence of shallow lattices, when the gap between sec- 
ond and third band is negligible. For larger s, such single 
particle excitations travel with smaller velocity, disperd 
more quickly, and can interfere with the propagation of 
the sound wavepackets. For this reason it is better to 
use slower perturbations T p 1 in order to restrict the 
dynamics to the first band. 



The shape of the wavepackets depends on the exci- 
tation producedure. However, the main features of their 
behaviour are quite general. In absence of the lattice, the 
observation of sound signals in a harmonically trapped 
condensate was achieved experimentally by employing 
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FIG. 1: GP simulation of wavepackets produced by a fast 
(T p = 1) perturbation of the type (j^Jl in presence of an optical 
lattice (s = 1, gn — 0.5Er). Both the first, the second and the 
third bands are excited leading to the formation of two pairs 
of wavepackets. For a slow perturbation (T p 2> 1) only the 
slow pair composed of phonons of the lowest band is created. 
Here, we plot the relative density n(x,t)/n(t,Q) as defined in 
Sect EH 

two different excitations methods ^]|: raising a poten- 
tial barrier in the center of the trap produces a density 
bump, which splits into two bright sound wavepackets; 
alternatively, removing a potential barrier from the cen- 
ter of the condensate gives rise to a dip in the density 
which splits into two grey sound wavepackets. 

The excitation method we adopt here is a combina- 
tion of these two: The initial condensate is in the ground 
state of a one-dimensional optical lattice. We then switch 
on and off a gaussian potential in the center and get two 
composed bright-grey sound signals propagating symmet- 
rically outwards. This procedure has the advantage that 
the ground states of the initial and final potential are 
identical. 

The perturbing potential has a spatial and temporal 
dependence 

Vp(x,t) = Vp x (x)V Pt (t), (3) 

where 

V Px (x) = bE R e^y[~x 2 /(wdf] , (4) 

!0, for t < 0, 

sin4 (llf 1 ) . for < t < T p h/E R (5) 
0, " for t > T p Ti/Er. 

The tunable dimensionless parameters are the width of 
the potential w, its height b and the time of perturbation 
T p . We impose the constraints 

w > 1 , (6) 

in order to address only the quasi-momenta in the central 
part of the Brillouin zone, and 

T p > 1 (7) 
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in order to excite the lowest Bogoliubov band only. Note 
that for typical densities, scattering lengths and lattice 
spacings, £ is of the order of d. Hence, w ^> 1 automat- 
ically implies wd S> £, which ensures that the produced 
excitations are phonons. 



The typical signals obtained in a GP simulation in the 
linear regime are shown in Fig. 1(2} , where in the upper 
planel we plot both the relative density n(x, t)/n(x, 0) 
and the convoluted signal. We denote by An and A</> 
the signal amplitudes of the density and relative phase 
wavepackets respectively. 



II. GROSS-PITAEVSKII (GP) AND DISCRETE 
NONLINEAR SCHRODINGER EQUATIONS 
(DNLS) 

We study the sound propagation in the presence of 
the lattice potential with the one-dimensional GP- 
equation 



ihtp 



2m 



V(x,t) + gnd\ip{x,t){' 



<fl(x,t), (8) 



describing a system uniform in the transverse directions 
y, z with n the 3D average density. The consequences 
of a transverse confinement, which can modify the de- 
gree of nonlinearity of the system and the speed of sound 
|28|. are discussed in Sect[V] The normalization of the 

wavefunction is f N ^ d ^ 2 \(p(x,t)\ 2 dx = N w , with N w 
the total number of lattice wells. The external poten- 
tial V(x, t) is given by the sum of the lattice potential 



opt 



sEji sin (nx/d) and the time-dependent pertur- 



bation Vp(x,t) written in Eq.(J3J) 

V(x,t) = V op t{x) + V P {x,t). 



(9) 



The ground state of the system in the presence of the op- 
tical potential is governed by two parameters: the optical 
lattice depth s and the interaction strength gn. 

Writing <p(x,t) = y/n(x,t)/n exp[iS(x, t)] , the GP- 
equation JSJl can be recast in the form of two coupled 
equations for the density and phase variables 



h(x) = 
S(x) = 



-d x 
1 



n(x) — d x S(x) 



Y{x,t) + gn(x)+ 



(10) 



■^V-n + ^S(x)f_ 



In the presence of the lattice, both the density and the 
phase are characterised by very strong modulations on 
the length scale of the lattice spacing d. In order to high- 
light the density variation corresponding to the sound 
wavepacket, it is in general convenient to plot the den- 
sity relative to the ground state density n(x,t)/n(x,0). 
Alternatively, in order to mimic the limited resolution 
of a detection system, one can perform a convolution 
over a few lattice wells of the GP results for n(x,t). We 
also consider the evolution of the relative phase between 
neighboring wells, defined as 



W1/2 = S(x =(£+ l)d) - S(x = id). 



(11) 
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FIG. 2: Typical result of a GP simulation for small perturba- 
tions (linear regime). Two sound packets move simmetrically 
outward at the Bogoliubov sound velocity. The relative den- 
sity n{x, i)/n(x, 0) (thin line), the convoluted signal (thick 
line) and the phase difference <^+i/2 are shown. The signal 
amplitudes of the density and relative phase are denoted by 
An and A<f) respectively. 



In the tight binding regime, where the condensate 
phase is approximately flat in each well, the GP dynam- 
ics can be simplified by using the discretised wavefunc- 
tion ipi = y/n e (t)exp[iS e (t)], where n e (t) = \ijje(t)\ 2 = 

fu-d/2 \^p{x,t)\ 2 dx and Si{t) = S(x = £d,t). At equi- 
librium ni = 1 and Si — 0. The time evolution of ip£ 
is given by the discrete nonlinear Schrodinger equation 
(DNLS) [Hll 

s 

ihipt - -- (^+i + + [Vt(t) + u\Mt)\ 2 ] Mt)- 



(12) 



Here, the external potential Vi(t) includes only the per- 
turbation Vp(ld,t), with Vp given by Q, since the pres- 
ence of the optical lattice potential is accounted for by 
the discretization of space. The two quantities 8 and U 
entering the DNLS describe, respectively, the tunneling 
coupling between two neighbouring wells (corresponding 
to half the height of the lowest Bloch band) and the on- 
site interaction U (corresponding to the inverse compress- 
ibility) . 

In terms of density and phase variables, the DNLS 
equation becomes 



J2 T^n t {t) ni ,{t) sm[St{t)-S t ,(t)], (13) 



t'=t±\ 
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S e = 



ln e ,(t) 
n e (t) 



cos[S e {t) - Se>(t)} 



1.5r 



where \ii = Uni + Vt. 

The results of the numerical integration of the DNLS 
equation l|12|) match those obtained by solving the GP- 
equation |JSJ in deep lattices. In order to compare the 
density n(x, t) calculated with the GP equation with 
the solution ne(t) of the DNLS equation, it is neces- 
sary to average n(x, t) over each lattice site. Instead 
the comparison for the evolution of the relative phase 
<Pe+i/2 = Sfr+i — Se is strighforward. 

The DNLS approach is not only convenient from the 
numerical point of view, but also allows a clearer and 
semi-analitycal understanding of the basic physics under- 
lying sound propagation and shock-waves formation in 
the tight binding regime. In particular the DNLS equa- 
tion written in the form (|13|) is particularly useful to get 
insight into the breakdown of the linear regime for deep 
lattices (see discussion in Sect II V|) . 

In the following, we will discuss the solution of the GP- 
equation (jHJ and DNLS equation (|12|) for various values 
of the potential height 6, the potential width w and the 
perturbation time T p . We will vary the parameters s and 
gn in the case of the GP equation and the parameters 
5 and U in the case of the DNLS to explore different 
regimes of lattice depths and interaction strengths. 



III. VELOCITY OF PROPAGATION OF 
WAVEPACKETS 
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FIG. 3: GP simulation for the relative density n(x, t)/n(x, 0) 
at t — 480h/En with gn — 0.5Er for different lattice depths: 
s = (a), s = 10 (b) and s = 20 (c). 




The first quantity we want to extract from our simula- 
tions is the velocity of propagation of sound wavepackets 
as a function of lattice depth. The result in the linear 
regime has already been predicted in |29j. Here we fo- 
cus on the case of large perturbations producing signals 
which can be observed experimentally. 



In Fig. |31 we plot some examples for the relative density 
at the final stage of our GP-simulation at lattice depths 
s = 0, 10, 20 with gn = 0.5-Er for a large perturbation. 
Those signals involve significant nonlinear effects, even in 
the uniform case (see FigOJa)). From these simulations, 
after performing the convolution of the signal, we obtain 
the sound velocity shown by the circles in Fig. The 
respective signal amplitudes An, together with the Bo- 
goliubov prediction (solid line) are also indicated in the 
figure. 

The first conclusion is that from our simulation we can 
extract the value of the Bogoliubov sound velocity with 
high accuracy. Moreover, we find that sound signals of 
measurable amplitude can be observed also in deep lat- 
tices where the sound velocity is considerably lower than 
in the uniform system. At this value of gn, relatively 
large signal amplitudes can be obtained up to s = 20. 



FIG. 4: Sound velocity as a function of lattice depth s. Bo- 
goliubov prediction (solid line) and results "measured" based 
on the simulation (circles) with respective signal amplitudes 
An for gn = 0.5Er. The signal amplitude An is defined as 
indicated in Fig. [5] 

Note that the signal amplitudes obtained from the simu- 
lation at s — 20, 30 correspond to their saturated values 
as discussed in detail further below. 



IV. LINEAR AND NONLINEAR REGIMES 

By keeping lattice depth and interaction fixed while 
increasing the strength of the external perturbation, the 
role of nonlincarities becomes more and more important. 
The system passes through three regimes: 

1. linear regime, where the Bogoliubov description 
holds and the variations of density and relative 
phases are small; 

2. shock wave regime, where density variations induce 
mode-coupling among Bogoliubov excitations, giv- 



5 



ing rise to the formation of shock waves. Depend- 
ing on the curvature of the Bogoliubov dispersion 
beyond the phononic regime, shock waves emerge 
in front of the wavepacket (uniform system, shallow 
lattice), both in the front and in the back (interme- 
diate lattice) or only in the back (deep lattice with 
8 <C U/2>). In the first two cases, the sound sig- 
nal is deformed and it is finally disperded. Instead, 
when shock waves form only in the back, the signal 
maintains a compact shape and propagates at the 
sound velocity predicted by Bogoliubov theory. 

3. Saturation regime, where the sound signal am- 
plitude saturates. The sound wavepacket leaves 
behind a wake of noise, and still propagates at 
the sound velocity predicted by Bogoliubov theory. 
This regime exists only in the presence of the lattice 
and provided that 8 -C U/3 to ensure that shock 
waves form in the back of the signal as explained 
in regime (2). For fixed gn or U, the condition 
8 <C U/3 can be ensured by making the lattice suf- 
ficiently deep. 

The signal amplitudes attainable in each regime and the 
perturbation parameters b, w, T p needed to reach a cer- 
tain regime depend on the lattice depth s and on the in- 
teraction strength gn, or equivalently, on 8 and U . In the 
presence of a lattice, a stronger perturbation is needed 
to obtain the same signal amplitude as in absence of lat- 
tice. This reflects the fact that the condensate in a lattice 
is less compressible. In the following subsections we are 
going to discuss more in detail the three regimes. 



1. Linear Regime 

In the linear regime, the signal moves without disper- 
sion and with sound velocity predicted by the Bogoliubov 
theory. The amplitudes An and A(f> are small and remain 
constant during the propagation. The ratio between the 
relative density and the relative phase amplitudes can be 
derived analytically. This can be done using the DNLS 
equations |(T3)l valid only for deep lattices, or the macro- 
scopic hydrodynamic formalism developed in |29| valid 
for all lattice depths and small quasimomentum. 

A solution is provided by the ansatz 

n(x,t) = I + An[f+(z - ct) + f-(x + ct)], (14) 
<f>(x,t) = -A<j>[f+(x - ct) - f-(x + ct)], (15) 

where fi represent either the DNLS density or the macro- 
scopic averaged density. The functions / + and /_ de- 
scribe respectively the wavepackets moving to the left 
and to the right, and An and Acj> are their amplitudes, 
as shown in Fig|3 We obtain 



An 



1 



Urn 



-A0 



DNLS 




(16) 



relation, we have replaced sin[Se(t) — S#(t)] with its ar- 
gument, which implies Se(t) — Sfj(t) <C tt. When this 
condition is not satisfied, non linear effects due to the 
optical lattice becomes important, and the linearised for- 
malism breaks down. 



2. Shock Wave Regime 

The peculiarity of this region is the formation of 
shock waves. In the uniform front wave emits 

shock waves in the forward direction (see FigJSfa)). The 
stronger is the external perturbation, the stronger is the 
deformation and spreading of the sound wavepacket due 
to the emission of shock waves. We smoothen out the 
large density fluctuations in the shock wave region with 
a gaussain convolution of the signal over few lattice sites, 
mimicking the limited resolution of a detection system. 
A measurement of the position of maximum or the min- 
imum of the resulting signal, propagating, respectively, 
faster and slower than sound, yields a velocity which de- 
viates from the Bogoliubov prediction. However, in spite 
of the deformation of the signal, it is possible to extract 
the Bogoliubov value of the sound velocity as shown in 
Fig0| by following the center-of-mass position of the sig- 

nal 113- 
In a shallow lattice, shock waves form in the front as 
in the uniform case, because the formation of a gap in 
the Bogoliubov spectrum does affect only a small range 
of quasi-momenta close to the Brillouin zone boundary. 
Hence, the mode-coupling among Bogoliubov excitations 
leads to the creation of excitations outside the phononic 
regime which travel at a speed larger than the sound 
velocity. 

On the other hand, in intermediate and deep lattices 
shock waves can also form behind the sound packet due to 
the behaviour of the lowest band Bogoliubov dispersion, 
which in the tight binding regime takes the form 0, |29| 



huj n 



1 25sin 2 



nq 
2^ 



2(5si 



nq 
2^ 



2U 



(17) 



where the last expression is written in terms of DNLS 
parameters. We notice that in order to get the DNLS 



For typical values of the density, the ratio 8/U lies be- 
tween zero and one. If the ratio 8/U is larger than 
1/3, the Bogoliubov dispersion has a positive curvature 
in a small range of quasimomenta and becomes nega- 
tive closer to the zone boundary. In this case the shock 
waves are produced both in the front and in the back 
of the sound wavepacket. The dispersion has a nega- 
tive curvature for all q as long as 8/U < 1/3. In deep 
lattices (where 8/U -C 1/3), wavepackets composed by 
quasi-momcnta outside the phononic regime will propa- 
gate much slower than sound and shock waves are there- 
fore created behind the sound wavepacket (Fig[3{b) and 
FigEJl. In our simulations, for 8/U = 10 -2 , we observe 
that the relative phase distribution is strongly deformed 
if <fie+i/2 ~ 7r /2, due to the non trivial sin-dependence 
of the current in the DNLS Ea.lfT3)l (see Fig. (upper 
two panels)). However this behaviour does not become 
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FIG. 5: DNLS simulation of the density n£(t) and relative 
phase 4>i+i/2 for a deep lattice (S/U = 10 -2 ) in the shock 
wave regime (T p b/w ~ 5.5). Upper panels: early stage (t = 
16h/En); lower panels: late stage (t — 200ft/ 'Er). 



FIG. 6: DNLS simulation of the density ne(t) and relative 
phase 4>e+i/2 for a deep lattice (S/U = 10 -2 ) in the saturation 
regime (T p b/w ~ 7). Upper panels: early stage (t = lGh/En); 
lower panels: late stage (t = 200H/ Er). 



critical up to the point where the relative phase of two 
neighboring lattice sites <j>e+i/2 = it, which defines the 
onset of regime (3). 



3. Saturation Regime in deep lattices 

If we further increase the strength of the external per- 
turbation, the relative phase 0^+i/2 at some site reaches tt 
(see FiglU upper two panels) and there the current starts 
flowing in the opposite direction (see Eq. (|13[) L As a con- 
sequence, a wake of noise is left behind the sound packets 
and we find a saturation in the amplitude of the propa- 
gating signal (see Figl^c) and FigEl lower two panels). 
The interesting feature is that the noise has zero average 
velocity, since the oscillations of population between dif- 
ferent wells are completely dephased. Provided that the 
lattice is deep enough to satisfy 5 <C U/3, shock waves 
form only in the back of the signal, as discussed above. 
Hence the noise and the shock waves never overtake the 
sound signal, which is always able to "escape" from them. 

We stress that this effect appears only in presence of a 
deep lattice. In fact, in the uniform case, even in presence 
of strong nonlinearities leading to a strong deformation 
of the signal, in the central region the system is always 
able to recover its ground state after the sound wave has 
passed by (see Fig|3](a)). 



To demonstrate the saturation of the signal, we follow 
the evolution of the system for a long time and look at the 
amplitudes An and A(j>. Both in the GPE (for relatively 
deep lattices) and in the DNLS simulations, we find an 
interesting scaling law that helps to distinguish between 
the three regimes mentioned above and makes evident 
the saturation of the signal amplitude. The scaling of 
the results is shown in Fig[7| The effect of the perturbing 
potential takes a universal form when: 

• the perturbation parameters are combined in the 
form Tpb/w which reflects the capability of the sys- 
tem to react to an external perturbation; 

• the relative density variation is rescaled as 
Any^U/S, while the amplitude of the relative phase 
signal A(/> does not need to be rescaled. 



The results summarized in Fig[7| show that for small 
parameter T p b/w, the perturbation produced in the sys- 
tem is small, and depends linearly on T p b/w. Increasing 
the perturbation parameter T p b/w, the signal amplitude 
An quickly saturates. The three regions indicated in 
Fig0 correspond to the three regimes (f-3) listed above. 

From our numerical results, we find that the propor- 
tionality between the relative density and the relative 
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FIG. 7: Results for the density and relative signal amplitudes 
An and A(f> obtained from DNLS simulations for varying per- 
turbation parameters T p , b, w (see Eq.©), tunnel coupling 8 
and on-site interaction U (see Ea, 1112(1 *1 with S <C 17/3. 



phase amplitudes written in Ea. (|16|) . holds also beyond 
the linear regime. In particular, we found a saturation 
of the relative phase amplitude at a value A<fi « 0.2ir, 
which goes with the saturation of the density amplitude 
An. This implies that for given interaction and lattice 
depth, the amplitude of the density signal is limited by 
approximately An « 0.2tt^S/U. 

Saturation does not occur in the uniform system or 
at low lattice depth. To demonstrate this, we plot in 
Fig. [SI the signal amplitude An measured after different 
propagation times for s = (dashed lines) and s = 15 
(solid lines) as a function of potential height b, with fixed 
potential width w and perturbation time T p . With s = 
15 the signal amplitude takes exactly the same values 
at t = 100, 200, 300, 400h/E R (the corresponding four 
lines in Fig. [^perfectly overlap!). In contrast, the signal 
amplitudes in the uniform system (s = 0, dashed lines) 
at different times (t = 100, 200Ti/Er) differ from each 
other: At earlier times the signal amplitude increases as 
a function of the strength of the external perturbation 
and does not saturate. The amplitude measured at later 
times coincides with the one measured earlier only for 
small potential heights 6, when the propagation dynamics 
is linear. When nonlinear effects becomes important, An 
decreases during the evolution as a consequence of the 
dispersion of the signal. 



s 
< 



0.2 



T b/w 
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FIG. 8: The signal amplitude An as a function of the per- 
turbation parameter bT p /w at gn = 0.5Er. Dashed lines: 
GP simulation for s = and gn = 0.5-Er at t = IOOTi/Er 
(upper dashed line) and t — 2QQTi/Er (lower dashed line). 
Solid lines: GP simulation for s = 15 and gn = 0.5Er at 
t = 100, 200, 300, 400ft/£ fl for s = 15. The corresponding 
four lines exactly overlap. 



V. ROLE OF HARMONIC TRAPPING 

Realistic experimental setups include an harmonic 
trapping potential in the parallel and transverse direc- 
tions of the lattice. 

In the longitudinal direction, the harmonic trap has a 
negligible effect on sound propagation if L S> wd, being 
L the size of the condensate and w the width of the per- 
turbing potential. This condition guarantees, first of all, 
that the sound wavepacket travels in a region of approx- 
imatively constant average density. Moreover it ensures 
that low energy discretised modes [l5l are not strongly 
excited. 

In the Thomas-Fermi regime, we find that the time 
required by the sound signal to travel along a distance 
of the order of the size of the condensate is ~ 2tt/u)d, 
with ujp — \J m/m*ui x the renormalised dipole frequency 

The transverse harmonic trapping leads to an inhomo- 
geneous radial density profile which can affect the sound 
velocity. Even in the uniform radial Thomas- 

Fermi density distribution ch ang es the sound velocity c = 
y/U/m to c = v / U/2m [H E3 Q E3 . The analogous 
modification of the velocity of sound to c = ^JU /2m* is 
expected in presence of the lattice provided the inverse 
compressibility U is linear in the density, as discussed 
in |28l l29l l30| . This approximation requires the inter- 
actions gn to be sufficiently small and has proven ap- 
plicable to describe the experimental result on collective 
excitations [34) and on the condensate size [3f| where 
gn « 0.2 ~ 0.5E R . 

Moreover the predicted effects require single band dy- 
namics. This is ensured if the radial trapping frequency is 
larger than the width of the first Bogoliubov band. In the 
tight binding regime, this condition reads hcj^ > 2^/U5. 
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VI. CONCLUSIONS 

In conclusion, wc find that there exists a broad range 
of optical potential depths where the change in sound ve- 
locity induced by the lattice should be measurable also 
in the non linear regime. For instance, a velocity of prop- 
agation c = 0.65c s=0 and 0.30c s= o can be measured at 
an optical lattice depth of s — 10 and 20 with a maximal 
density variation An max = 0.13 and 0.05 respectively. 
This maximal attainable signal amplitude decreases very 
strongly with the optical lattice depth, making the signal 
in practise observable only up to a certain lattice depth 
which depends on the interaction strength. 

We have shown that the presence of a deep lattice has a 
dramatic effect on the propagation of sound signals in the 
nonlinear regime: Contrary to the uniform system, shock 
waves propagate slower than sound waves, due to the 



negative curvature of the dispersion relation in the lowest 
Bogoliubov band. Moreover, nonlincarities can play a 
role also at very small density variations and induce a 
saturation of the sound signal, which goes along with 
dephased currents left beyind the signal. This effect has 
no analogue in the uniform case. One should also keep in 
mind, that the saturation effect is clearly evident only if 
the shock waves move with a velocity much lower than the 
sound velocity, which requires S/U <C 1/3. Obtaining a 
clear saturation effect (small S/U) associated with a large 
signal amplitude An (large S/U) at fixed lattice depth 
requires a compromise in the choice of lattice depth s 
and interaction strength gn. 
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